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^~, — ■ ^ Abstract. Starting from Newton's equations of motion, we derive a dynamical 

Q ■ density functional theory (DDFT) applicable to atomic liquids. The theory has the 

"^ , feature that it requires as input the Helmholtz free energy functional from equilibrium 

density functional theory. This means that, given a reliable equilibrium free energy 
functional, the correct equilibrium fluid density profile is guaranteed. We show that 
when the isothermal compressibility is small, the DDFT generates the correct value 
for the speed of sound in a dense liquid. We also interpret the theory as a dynamical 
equation for a coarse grained fluid density and show that the theory can be used 
O I (making further approximations) to derive the standard mode coupling theory that 

is used to describe the glass transition. The present theory should provide a useful 
starting point for describing the dynamics of inhomogeneous atomic fluids. 
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'nJ" ■ 1. Introduction 

o 

O '. Classical density functional theory (DFT) jTJ |2] allows one to determine the equilibrium 

ensemble average one-body density profile of a fluid in the presence of an arbitrary 
external field \^^^*(r). As DFT is based upon a minimisation principle for a free energy 
density functional [HE], one can use DFT to calculate thermodynamic quantities such as 



^ [ surface tensions. In practice, DFT is not exact because the exact free energy functional 

is generally unknown. Nonetheless, the theory has proved to be very powerful and 



there exist in the literature a large number of rather accurate approximations for the 
'k/< I Helmholtz free energy functional [2j and equilibrium DFT is a well developed theoretical 

H I framework for describing static fluid phenomena jB]. Given this body of knowledge, 

constructing a dynamical density functional theory (DDFT), that requires as input the 
same Helmholtz free energy functional as the equilibrium DFT, for the dynamics of 
the one-body density profile of an inhomogeneous fluid that is out of equilibrium, is a 
very appealing idea. In recent work |3I Ej , Marconi and Tarazona derived an equation of 
motion for the one body density profile of a fluid of particles with stochastic equations of 
motion (i.e. for colloidal fluids). This theory has been applied to a number of problems 
[HI 13 El El Hm HH 1121 Uni mi Uni and has proved to be rather accurate where comparison 
with the results of Brownian dynamics simulations has been made. The question that 
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then naturally arises is, can one construct a DDFT for a fluid of particles with motion 
governed by Newton's equations? In a recent paper, Chan and Finken ^Hl proved that 
the non-equilibrium one body density p(r, t) is uniquely determined for a given time 
dependent external held l^'^^*(r, t), where t is the time. They showed that in general 
the DDFT requires a non-equilibrium free energy functional that is not the same as the 
free energy functional entering equilibrium DFT and that the (approximate) theory of 
Marconi and Tarazona [HEj can be obtained by taking the adiabatic limit of the 'exact' 
DDFT jin]. Given that very little is known about the non-equilibrium free energy 
functional, this makes practical calculations of the fluid dynamics in the Chan-Finken 
framework somewhat difficult. 

In the present paper, we derive an approximate DDFT for a fluid of particles 
described by Newton's equations of motion (i.e. an atomic fluid, as opposed to a colloidal 
fluid). Our DDFT requires as input the equilibrium free energy functional and therefore 
makes practical calculation of the fluid dynamics much more tractable. The theory is 
an approximate one, and in the form derived here is most relevant to dense fluids. 
On applying the theory, in the limit that the ratio of the heat capacities CpjCy -^ 1 
(i.e. when the isothermal compressibility is small), we flnd that the DDFT reproduces 
the correct value for the speed of sound waves in a dense liquid, and predicts the length 
scale over which sound waves are attenuated. We also interpret the theory as one for 
a coarse grained one-body density proflle, p(r, t), and from these equations we derive 
(making a number of further approximations) the mode coupling theory (MCT) that 
is used to describe the glass transition [HI Ej . We beheve that the present theory will 
provide a useful starting point for determining the dynamics of the one body density 
profile of an atomic fluid. 

This paper is laid out as follows: In Sec. |2l starting from Newton's equations of 
motion, we derive equations of motion for the ensemble average fluid one body density 
profile p(r, t). In Sec. 01 we introduce a number of approximations in order to close the 
equations derived in Sec. El This closure yields us the DDFT. In Sec. 01 we apply the 
DDFT in two particular situations: (i) sound waves propagating in a fluid and (ii) to 
dynamic two point correlation functions, where we derive the MCT. Finally, in Sec. El 
we draw some conclusions. 

2. Density equations of motion 

The (classical) density operator for a system of N identical particles is /3(r, t) = 
^^(5(r — rj(t)), where r^ = {rj(t)} is the set of position coordinates of the particles and 
the index i labels the different particles. In Fourier space, the density operator is given 
by: 

Pk(t)= /dre^''-p(r,t) = 5^e^''-^ (1) 
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Newton's equations of motion for these particles are 

m OYi 
where m is the mass of a particle and 

y(r^,t) =Y,V^^\v,,t) + \Y,v{v,{t) - r,(t)) (3) 

is the potential energy of the system. We assume this is composed of one body 
contributions due to the external potential l^^^'*(r, t) and a sum of pair interactions 
between the particles; v{r) is the pair potential. Following Zaccarelli et al [IHj, we take 
two derivatives of Eq. (^ to obtain: 

Pk(t) = - 5^(k ■ r.(t))^e*^-^» + J2'^^- r^W)^'''"'- (4) 

i i 

We can express the potential energy as: 



«j 



in an obvious notation. Differentiation gives: 

Combining Eqs. Q, Q and © we obtain the following equation |18j : 
M) = - E(k ■ '.('))'<=*"■" - ;^ E / T^S*!'' ■ k')e-'-''-'-'')e--- 

-iHJ (^^-"W(k ■ k')e--'-e--. (7) 

Using Eq. (P) this can be rewritten as: 
hit) = - $^(k ■ r.(t))V^-r« - 1 1 7^^k'(k ■ k')pk-k'(t)pk'(t) 

dk' 



my (2^)3^i^"W(k-k')Pk-k'(t). (8) 

This equation is simply a generalisation to the case with an external field of Eq. (8) in 
Ref. {18 . We transform Eq. (jH)) back into real space and obtain the following: 

/^(r, ^) = - E / (|^(k ■ r,(t)) V'^-(r-r) + ;^^ ■ / ^''^^^' ^)'^^''' ^^^'^^^ " ^'^ 

+lv.[p(r,t)VV^^'*(r,t)]. (9) 

m 

Our interest is in ensemble average quantities. We therefore average in Eq. over 

the ensemble of initial momenta and positions to obtain an equation of motion for the 
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ensemble average one body density, p{r,t) = (p(r, t)): 

p(r, t) = -(^ I ^(k • r,(t)) ^'■(r-r)) + ^ V ■ I dr'p(') (r, r', t)Vv{r - r') 

+ -V-[p{r,t)W'\r,t)], (10) 

m 

where (■ ■ ■) denotes the ensemble average and where p^'^\r,r',t) = (p(r, t)p(r', t)) 

is the ensemble average, time dependent, two body distribution function. We now 

consider the first term on the right hand side of Eq. (jlOj] . This contains the expression 

k ■ rj(t) = kr-i^parit), where ri^parif) is the component of the velocity of the i^^ particle in 

the direction of k. We write rlparit) = {'^f^par) + (^iif)^ where aiit) is just the fluctuation 

of rf^parif) about its mean equihbrium value. Note also that the equilibrium quantity 

(r| ) = kBT/m, from equipartition. We may therefore re-write the first term on the 

right hand side of Eq. ()1()|1 as follows: 



E / (^"^ ■ *. WJ^-^"-"--") = / (|^*'<E ".(*)e-->e- 



ikr 



ksT / dk , 2/V^ ikri\„-Jk-r 



= S(r,t)-^VV(r,t). (11) 

In order to obtain the second line we have used Eq. (Q) and defined 

B{r,t) ^ I ^,k^J2^,{t)e^'"^^)e-^^-\ (12) 

Substituting Eq. (fTT|) into Eq. (fTUI) we obtain 

p(r, t) + B{r, t) = ^VV(r, t) + -V ■ f dr'p^^\r, r', t)Vv{r - r') 
m m J 

+-V- [p(r,t)Vr^"*(r,t)l. (13) 

m 



This equation, taken together with Eq. ()12|1 is still formally exact. We note that the 
right hand side of this equation is simply the derivative of the first equation of the YBG 
hierarchy j3] . This means that at equilibrium the right hand side of Eq. fll3|) is equal to 
zero and therefore at equilibrium B{r,t) = 0. In the Appendix we present an alternative 
derivation of Eq. ()13p. starting from the first equation of the BBGKY hierarchy |2]. The 
derivation gives us an alternative expression for B{r, t), that may provide further insight 
into the nature of this term. 

3. The DDFT 

Eq. (fT^ is the starting point for our derivation of the DDFT. We make two 
approximations in Eq. (fT!?|) in order to render a closed and tractable DDFT. Our first 
assumption is the same as that used to derive the the colloidal DDFT [U [71 |H] : we 
assume that the two particle correlations in the non-equilibrium fluid are the same as 



Dynamical density functional theory for dense atomic liquids 5 

in an equilibrium fluid with the same one body density profile, i.e. we assume that we 
can apply the following sum rule, which is exact for an equilibrium fluid [7j: 



- kBTp{v)Vc^^\r)= f dr'p^^\r,r')\/v{r,r'), 



(14) 



where c^^\r) is the one body direct correlation function and is equal to the functional 
derivative of the excess (over ideal) part of the Helmholtz free energy functional ^ |2] : 

,a,,r) . -/-^. (15) 

We must also make an approximation for the term B{r,t) in Eq. |T!?|) . Since at 
equilibrium B{r,t) = 0, the simplest approximation we can make is to assume 
B{r,t) oc p(r, t). Making these two approximations we obtain our main result: 



d'pir,t) ^ ^ gp(r,t) _ 1^ 
dt"^ dt m 



p(r,t)V 



SF[p{r,t)] 
6p{r,t) 



(16) 



where u is an undetermined (collision) frequency. In general, we expect i/ to be a (non- 
local) functional of the density p{r,t). We show below that when density modulations 
are small, then the collision frequency u = ksT /niD, where D is the self diffusion 
coefficient. F is the Helmholtz free energy functional: 



F[p(r, t)] = FuW, t)] + Fe,[p(r, t)] + j drV;,i(r, t)p(r, t). 



(17) 



where 



F,,[p(r,t)] = ksT J drp{r,t)[\n{p{r,t)A^) - 1] (18) 

is the ideal gas contribution to the free energy, and where A is the thermal de Broglie 
wavelength. We must emphasise that strictly speaking, F is defined as a functional 
of the equilibrium ensemble average density density profile p(r). We are effectively 
assuming that we may use F as an approximation for the unknown dynamical free energy 
functional. It is important to note that at equilibrium, given a reliable approximation for 
F, Eq. (fTH|l guarantees the exact equilibrium fiuid density profile [Zj. For the case when 
the potential energy l^(r^,t) also contains three-body and higher body interactions, as 
we show in the Appendix, one still obtains Eq. (fTT)|l . making approximations equivalent 
to those made above. The argument used is somewhat similar to the approach used 
in Ref. [7] to generalise the DDFT of Marconi and Tarazona [H |3] to the case with 
multi-body interactions. 

In order to show that when density modulations are small the coefficient u = 
ksT /niD, we consider a fiuid in which a (small) number, n, of the particles are labelled 
in such a way that we can think of them as a different species of particles, without 
changing any of the interaction potentials between the particles. If at time t = the 
n labelled particles are confined to a small volume near the origin (r = 0), then in the 
limit t ^ oo the density profile of these particles will be pi{r,t) = nGs{r,t), where / 
denotes the labelled particles and Gs{r,t) is the van Hove self correlation function [Hj. 
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In this limit the one body density of the remaining particles p(r, t) — Pb, where pb is the 
bulk fluid density. It is known that in the limit t ^ oo that |3j 

Using this result, together with Eq. (fT^ generalised to a two component fluid, we 
obtain the following after making a Taylor expansion of the excess Helmholtz free energy 
functional: 

ot^ ot m 

where ■ ■ ■ denotes terms of order pf, pi{p — pb) and higher order terms. Furthermore, 
in the limit t ^ oo, p^ ^ p^, so we may neglect the pi term. Substituting Eq. (fT^ into 
Eq. ()2()|1 we obtain the result that u = ksT /mD. 

At this point, we note that in cases where the collision frequency v is large, then 
up ^ p and we may, in certain circumstances, approximate Eq. p6|) by 



dt m 



p(r,t)V.^^[^('^'^)] 



(21) 



6p{r,t) 

This is equivalent to the DDFT of Marconi and Tarazona jH |3] , which was derived for 
colloidal fluids (fluids of Brownian particles) - see also the discussion in Ref . ^\ . 

4. Applications of the DDFT 

Having argued for considering Eq. (fTBj) as an approximation to the equation governing 
the dynamics of the particles, we now examine some of the consequences. The flrst case 
we consider is that of sound waves propagating in a dense fluid. 

4-.1. Sound waves in a dense liquid 

We Taylor expand the excess Helmholtz free energy functional in terms of density 
modulations Sp{r, t) about a uniform bulk density pb (i.e. we write p(r, t) = pb + Sp{r, t)). 
Omitting terms beyond 0{Sp'^) we obtain J21IZ1: 

Fe,[p(r,t)] = Fe,[pb] -c«(oo) fdr6p{r,t) 

-^ I dr I dr'5p(r, t)Sp{r', t)c<^'^ (r - r') , (22) 

where c^'^\r) is the Ornstein-Zernike pair direct correlation function. Combining 
Eq. dig) with Eq. ^ we obtain ^: 

-KTo + ^ ^: = V5p(r,t) V / dr5p(r,t)Vc^ ^(r-r). (23) 

ot^ ot m m J 

Sound waves correspond to the following fluid one body density proflle: 

p(r,t) = pb + 5p{Y,t) 

= Pb + eexp[2(q-r-u;t)], (24) 
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where e is the amphtude, which is assumed to be small compared to the bulk fluid 
density p^ and |q| = k + i\~^ . The speed of sound in the liquid, Cg = uj/k, and A 
is the distance over which the amplitude of the sound is attenuated by a factor 1/e. 
Substituting Eq. ^ into Eq. (J2S1) we obtain: 

i-u^ - ^uu) 6p = (-'-^^ + fPl!t^lM] sp, (25) 

\ m m J 

where c(g) is the Fourier transform of c^'^\r). Equating the real and imaginary parts on 
the left and right hand sides of Eq. (j25j) we obtain the dispersion relation: 

and the following expression for the attenuation length: 

A(*) = |. (27) 

where we have used the fact that the fluid static structure factor S{k) = (1 — pbc{k))^^ 
and we have assumed that k ^ A^^. 

Sound waves correspond to the long wavelength limit. In this limit, S{k) ~ 5(0) = 
pijksTxT-, i-e. proportional to the isothermal compressibility xr, yielding 

cl = -^. (28) 

PhmxT 

The adiabatic speed of sound in liquids is known to be [Sj 

cl = ^— , (29) 

PbinxT 

where 7 = Cp/Cy, is the ratio of the heat capacities. The reason that our expression 

(OHj) for Cg differs from the true value, given by Eq. (j^ . can be traced back to the 

approximation B{r,t) ^ z/p(r,t), used to obtain Eq. (fT?)|l from Eq. pH|l . In making 

this approximation we are effectively neglecting temperature fluctuations. However, the 

density fluctuations corresponding to sound waves occur on time scales too fast to reach 

isothermal conditions; sound wave fluctuations are an adiabatic processes. Using the 

fact that Cp = Cy + TxT^v/Pb [3 , where Pv is the thermal pressure coefficient, we see 

that it is only at low temperatures and high densities, when xt is small, that Cp ^ Cy, 

and Eq. ()28|) yields the correct value for Cg, in agreement with Eq. ((221) • 

4.2. DDFT and MCT 

Using DDFT (in various formulations) to describe the slow dynamics of fluids at state 
points near to the glass transition has been suggested by Kawasaki and coworkers 
Pm EB 1221 123 121] ■ in this section we show how the standard mode coupling theory 
(MCT) can be obtained from a temporally coarse grained DDFT, as opposed to the 
spatially coarse grained DDFT, described in Ref. [20] • Our approach, in principle, 
allows us to include the effects of external potentials in the MCT. The arguments in 
this section are to some extent equivalent to those made in Ref. ^18., . 



-p''' 
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Eq. (fTH|). as derived above, is an equation for the ensemble average fluid one body 
density profile, p{r,t). However, we may also interpret it as an equation for the time 
evolution of a coarse grained density profile, p(r, t). For example, we define a temporally 
coarse grained density jB]: 

/oo 
dt'Kit-t')pir,t'), (30) 

-oo 

where K{t) is a normalised function of finite width in time. Multiplying Eq. Q by 
K{t — t') and then integrating over t, we obtain an equation of the same form as Eq. (|T0|) 
for the coarse grained density p(r, t). This equation involves the coarse grained two body 
density distribution function: 

/oo 
<lt'K{t~t')p{r,t')p{v',t'). (31) 

-oo 

If we assume further that we have coarse grained sufficiently that we can approximate 
the term involving p^^)(r,r', t) using the sum rule ()14|). and that the coarse grained 
equivalent of -B(r, t), i? ~ up, then we obtain Eq. (fT6|l as our equation of motion for the 
coarse grained density p(r, t). We are effectively assuming that we can approximate the 
coarse grained free energy functional by the Helmholtz free energy functional entering 
equilibrium DFT - i.e. the free energy F in Eq. ()17j) is now a functional of the coarse 
grained density p{jc,t). 

We now proceed to obtain the MCT, using a particular approximation for F, namely 
the Taylor expansion in Eq. ()22p . We will also assume for the time being that the 
external potential V^xt = 0. What follows is to some extent the argument of Kawasaki 
in Ref. [22] • Using Eq. (jSl in Eq. ^ we find (Zj: 

at^ ot m 



+^— V- 
m 



(Pb + 6p{r,t)) Jdr'6p{r',t)Vc^'\\r - r'|; p,)] = 0. (32) 



Fourier transforming Eq. (J32|l we obtain [2j: 

.2 „. l2 



k Ohk 11/" 

Pk(t) + z^Pk(t) = --—Pk{t) + ^c^p^it) + -— 7T^^ / dk'k ■ k'pk'(i)ck'Pk-k'(i), (33) 
pm pm pm [2tt)'^ J 

where Ck is the Fourier transform of the pair direct correlation function c*^^-* (r) , obtained 
by taking two functional derivatives of F[p{r, t)], and pk(^) denotes the Fourier transform 
of Sp{r,t). We rewrite Eq. (jH^ . in the form: 

Pk(t) + i^pUt) + fikPk(t) = Mt), (34) 

where ^2^ = fc^//5mS'k (we have defined the static structure factor ^k = [1 — p^Ck]^^) 
and where we have defined 

^k(i) = TT^^ /"dk'k • k'pk'(t)ck'Pk-k'(i^)- (35) 

Eq. (jM)) is suggestive, since it is known that one can always write formally the equation 
of motion for a dynamical variable such as Pk(i) in the following form |2SII2E]: 

Pk(t) + vp^,{t) + nlp^{t) = - f dt'mk(t')Pk(t - t') + i?k(t), (36) 

Jo 
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where Rk{t) is a noise term and the memory function is given formally by: 

,,. (-Rk(i)fi-k(O)) ,,_, 

where (■ ■ ■) denotes the average over the ensemble of initial configurations of the 
density. Since pk(^) = "^k ■ jk(^) = ikj^it), where j^if) is the longitudinal current, 
the denominator in (J37j) is simply equal to k'^NkBT/m. One finds that p_k(0) and 
R]i(t) are statistically independent |25j, so on multiplying Eq. (jHB|) by p-k(O), and then 
averaging, one obtains the following dynamical equation: 

Mt) + i^Mt) + ^iMt) = - [ dt'muit')Mt - t'), (38) 

where the normalised correlator 

m:> ^ p§^. (39) 

(Pk(O)p-k(O)) 
Eq. ()38|) is formally exact. However, one does not have a tractable form for the memory 
function m\^{t'). It is at this point that we return to the DDFT Eqs. (j34p and (j35|) and 
we make the following assumption j25j : 

(i?k(t)i?-k(0)) ~ (^k(t)/2-k(0)), (40) 

i.e. we will assume that the time correlations in -Rk(i^) and R]<i{t) are the same |25j . 
We also assume the usual MCT factorisation of four-point correlation functions into 
products of two-point correlation functions [SI EI: (Pk'(^)pk-k'(^)Pk"(0)p-k-k'(0)) — 
(Pk'(^)Pk"(0))(pk^k'(i)p-k-k'(0)) + (Pk'(^)P-k-k'(0))(Pk-k'(t)pk"(0)) and we obtain the 
standard MCT expression for the memory function [T7tl25]: 

^k(t) = ^!l''^!!\ I dk' [k ■ k'ck' + k • (k - k')ck-kf 5k'5k-k'0k'(i)0k-k'(i). (41) 

In this analysis we have derived equations for determining the correlator (jH^ . where 
Pk(^) is defined as the Fourier transform of the coarse grained density. However, 
the standard MCT is a theory for the correlator (jH^ . with Pk(t) given by Eq. ((T]). 
Strictly speaking, we have derived a theory with exactly the same structure, but for a 
different correlator. However, given that much of the physics of dense fluids near to the 
glass transition concerns the behaviour of collective density fluctuations, it is not too 
surprising that a description in terms of a coarse grained density is not that different 
from a description in terms of the density field itself. 

Returning to the case when there is an external field, we find there is the following 
additional term on the right hand side of Eq. ()34j) : 

Ek(t)=«k- /"drexp(zk-r)p(r,t)VV;xi(r,t). (42) 

This results in additional terms in Eq. (jHH|) as well as modifying the memory function 
mk(t). Consider an external potential Kxt(i') ~ exp(iq ■ r), i.e. a periodic field such as 
that generated by a laser with wavelength 27r/g. In this case we find 

Ek(t)~-ek-qpk+q(t), (43) 
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where e is proportional to the aniphtude of the external field. This term in Eq. (|34j) 
couples density fluctuations with wavelength k to fluctuations with wavelength k + q. 
Note also that the coupling is strongest when k is parallel to q. Including such a 
contribution in Eq. (jH^ . multiplying through by P-k(O) and then averaging would result 
in an equation relating the dynamics of 0k(^) to the correlator (pk+q(^)p-k(0)). For an 
equilibrium homogeneous fluid this correlator would be zero for q 7^ 0. However, for 
the inhomogeneous fluid this is not the case. This makes including the influence of an 
external fleld in the MCT a difficult problem! 

5. Conclusions 

In this paper we have presented a DDFT for systems of particles whose dynamics is 
governed by Newton's equations. In Sec. 121 we derived Eq. (|T3|l . an exact equation, 
which is the starting point for our DDFT. In Sec. El we suggested two approximations 
in order to close Eq. (fT^ . yielding the DDFT, Eq. ()16p. Then, in Sec. IHwe analyzed 
the consequences of assuming Eq. p(j|l . for the description of sound waves and for 
describing two point correlation functions. We have shown that the DDFT incorporates 
the physics described by the MCT, and the present analysis adds to the work of Kawasaki 
and coworkers [201 |2I1 122 I2B1 1211 1211 i^ suggesting that the DDFT provides a useful 
framework for studying the glass transition. 

We believe that Eq. (fT!^ needs further study. In particular, we believe that there 
may be better ways to approximate the term -B(r, t), given by Eqs. ()12p or ()56p . We made 
the very simple approximation -B(r, t) ~ z^p(r, t). However, when the modulations in the 
density profile become large we expect this approximation to break down. Whether this 
can be remedied by replacing the coefficient z/ by a function(al) of the density remains 
to be investigated. Having said this, the simple approximation suggested in Sec. El for 
the term B{r,t) in Eq. (fT^ . results in a dynamical theory that proves to be rather 
interesting, given that it can be used as a starting point for deriving MCT. 

Very recently, a DDFT (modified Cahn-Hilliard theory) of a similar form to that in 
Eq. (fTHjl was suggested in Ref. J2ZI. These authors considered a system of particles with 
stochastic equations of motion with correlations in the stochastic noise. If one considers 
a system of Brownian particles with Langevin equations of motion of the form 

r,(t) = -rV,y(r^, t) + rX,(t), (44) 

where F"^ is a friction constant characterising the one-body drag of the solvent on 
the colloidal particles and Xj(t) = {C,i (t) , ^f (t) , ^f (t)) is a white noise term with the 
property (C(^)) = and {^f{t)^j{t')) = 2kBT6ij6'"'6{t - t'), (where a,u = x,y,z label 
the Cartesian coordinates), then one can derive a DDFT of the following form 01|7|: 



dt 



p(r,t)V.^^t^(^'^)] 



(45) 



6p{r, t) 

However, if the noise term has memory - i.e. if the noise has the property (^f (t)) = 
and (C"(^)^j (^')> = 2kBT5ij5''''W{t - t'), where W{t - t') = exp(-|t - t'|/7)/7, then 
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arguing along the lines suggested in Ref. j27j, one arrives at the following (approximate) 
DDFT equation: 



^ dt^ dt 



P(M)V.^^[^^^'^)] 



(46) 



6p{r,t) 

The authors of Ref. J2ZI applied a linearised form of Eq. (jlBj) . together with a Ginzburg- 
Landau free energy functional for F to study spinodal decomposition. These authors 
were motivated by the need for a causality constraint in the description of phase 
transitions expected to occur in excited states of matter produced in heavy ion collisions 
|27j . They found that the memory term can have a significant effect in delaying the 
growth of density fluctuation in the initial stages of spinodal decomposition. Eq. (jlBj) 
has the same formal structure as Eq. ()l(j|l . We believe that it is worthwhile to pursue 
further investigations of these equations. In this sense, the present paper is a starting 
point for further developing a DDFT framework for tackling the microscopic dynamics 
of fluids. 
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Appendix 

Here we present an alternative derivation of Eq. ()13|). The state of a system of 
A^ identical particles is specified by the set of particle position coordinates r^ = 
{ri, r2, ■ ■ ■ , tn} and momenta p^ = {pi, p2, • ■ ■ , Pn}- There exists a phase space 
probability density function, f^^\r^, p^, t), which gives the probability that the system 
is in a particular configuration (r^, p^) at time t. The time evolution of /*^^^ is governed 
by the Liouville equation j2j: 

^ = m,/(^)}, (47) 

where {■, ■} denotes a Poisson bracket, and the Hamiltonian 

1 ^ 
7^(r^,p^t) = — 5^|p,p + r(r^t), (48) 

where V{r^,t), the potential energy of the system, is given by Eq. Q. Integrating 
over the Liouville equation, we obtain equations of motion for reduced phase space 
distribution functions /^"^ where n < N. One finds that in order to describe the 
dynamics of /'■"^ one must know member /("+^) of the BBGKY hierarchy j3]. The first 
member is the following: 



(l + |-^" + ^'-^-)^'"''-"P"') 
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= - j dv2 j dp2Fi2 ■ Vp J^'Hri, Pi, r2, P2, t), (49) 

where Xi = — Vri^'^'^*(ri, t) is the external force each particle 1 and F12 = — Vnf (ri — 
r2) is the pair force on particle 1 due to particle 2. If we integrate Eq. ^^ with respect 
to the momentum pi, then we obtain the continuity equation |28j : 

^ + V„.j = 0, (50) 

where 

p(ri,t) = y"dpi/«(ri,pi,t) (51) 

is the average one-body density and 

j(ri,t)= /dpiPi/«(r„pi,t) (52) 

J m 

is the current. If we now multiply Eq. ()49|1 through by Pi/m and then integrate with 
respect to the momentum pi, then we obtain the following momentum balance equation 

m- 

fi'j(ri,t) 



dt 



where 



+ Vr,- /dpiPi^/«(ri,pi,t)-lp(ri,t)Xi 

-- I Av2F^2p'^^\v^,V2,t) = Q, (53) 

m J 

p(2) (n, r2, t) = j dpi j dp2/(') (ri, pi, r2, P2, t) (54) 

is the two body density distribution function. Note that piPi is a tensor product 
(dyadic) J2H] • Taking a time derivative of Eq. (j^UI) and then using Eq. (j3^ to eliminate 
the term involving dj/dt we obtain the following equation: 

<9V(ri,t) ^ ^ /", PiPi.(i)/ .^ 

--Vr, ■ [ dr2FuP^'\r,,r2,t) - ^Vr, ■ [p(ri,t)Xi]. (55) 

m J m 

At equilibrium, from equipartition we find that the integral / dpi(pipi)/'^^^ = 
mkBTp{Ti,t)l, where 1 denotes the 3x3 unit matrix. Using this result in Eq. (jKKjl . we 
obtain Eq. (|T^. together with the following expression for B{r,t): 



Bir,t) = -V ■ V ■ /dp f^ - ^1) f''\r,p,t). 



(56) 



It is clear from this expression, that at equilibrium B{r,t) = 0. This analysis shows 
that B{r, t) is essentially the gradient of the off-diagonal elements of the kinetic part of 
the pressure tensor 
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We also note here, that if the potential energy V{y^ ,t) contains not just one and 
two-body contributions, as given by Eq. ©, but also contains three-body and higher 
body contributions: 

r(r^,t) = 5^ V-*(r„t) + ^ J]^2(r.(t),r,(t)) + ^ J]t;3(r,(t), r,(t), r,,(t)) + ■ ■ -, (57) 






then Eq. (j55|) becomes: 
9V(ri,t) 



9t2 



J 771^ 771 



--Vri- / dr2Fi2p('Hri,r2,t) 



-^Vr, ■ / dr2 y dr3Fi23P^'^ (ri, r2, r,3 , t) + ■ • • , (58) 

where F123 = — Vrif3(ri, r2, r3) is the three-body force on particle 1 due to particles 2 
and 3 and where 

p(=^)(ri,r2,r3,t) = / dpi / dps / dp3/(^Hri, pi, r2, P2, Ps, t) (59) 

is the three body density distribution function. In the case of many body interactions, 
the sum rule (fT^ becomes 

- fcBTp(ri)Vc«(ri) = j dv2P^^\vi,V2)VrM^i.^2) 

+ / dr2 / dr3p*^^)(ri,r2,r3)Vri^^3(ri,r2,r3) H 

= J] / dr2 ... / dr„p(")(r")Vr,t;„(r") (60) 

n=2 -^ ^ 

Applying this sum rule as an approximation for terms on the right hand side of Eq. (|58|) . 
together with the approximation B{r,t) ~ up{r,t), we obtain Eq. p6|) . 
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